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We study, using general scaling arguments and mean-field type calculations, the behavior of the critical 
Casimir force and its interplay with the van der Waals force acting between two parallel slabs separated at a 
distance L from each other, confining some fluctuating fluid medium, say a non-polar one-component fluid or 
a binary liquid mixture. The surfaces of the slabs are coated by thin layers exerting strong preference to the 
liquid phase of the fluid, or one of the components of the mixture, modeled by strong adsorbing local sur¬ 
face potentials ensuring the so-called (+,+) boundary conditions. The slabs, on the other hand, influence the 
fluid by long-range competing dispersion potentials, which represent irrelevant interactions in renormalization 
group sense. Under such conditions one usually expects attractive Casimir force governed by universal scaling 
function, pertinent to the extraordinary surface universality class of Ising type systems, to which the dispersion 
potentials provide only corrections to scaling. We demonstrate, however, that below a given threshold thickness 
of the system L cr j t for a suitable set of slabs-fluid and fluid-fluid coupling parameters the competition between 
the effects due to the coatings and the slabs can result in sign change of the Casimir force acting between the 
surfaces confining the fluid when one changes the temperature T, the chemical potential of the fluid p, or L. 
The last implies that by choosing specific materials for the slabs, coatings and the fluid for L < L cr i, one can 
realize repulsive Casimir force with non-universal behavior which, upon increasing L, gradually turns into an 
attractive one described by an universal scaling function, depending only on the relevant scaling fields related 
to the temperature and the excess chemical potential, for L i C rit- We presented arguments and relevant data 
for specific substances in support of the experimental feasibility of the predicted behavior of the force. It can 
be of interest, e.g., for designing nano-devices and for governing behavior of objects, say colloidal particles, at 
small distances. We have formulated the corresponding criterion for determination of T CI jt- The universality is 
regained for L 3> L cr j t . We have also shown that for systems with L < L cl -; t the capillary condensation phase 
diagram suffers modifications which one does not observe in systems with purely short-ranged interactions. 

PACS numbers: 64.60.-i, 64.60.Fr, 75.40.-s 


I. INTRODUCTION 

When a fluctuating field is confined by material bodies, ef¬ 
fective forces arise on them. This is due to the fact that the 
bodies impose boundary conditions on the medium, depend¬ 
ing on their geometry, mutual position and material properties, 
which leads to a modification of the allowed fluctuations in the 
medium. The last leads to a dependence of the ground state, 
or the thermodynamic potential of the system (say the free 
energy) on the geometry of the system and on the distances 
between its (macroscopic) components. In order to change 
these distances one has to apply a force that depends on the 
induced change of the allowed fluctuations. If the fluctua¬ 
tions are long-ranged the corresponding forces are also long- 
ranged. The existence of such long-ranged fluctuation medi¬ 
ated forces is called the Casimir effect and the corresponding 
forces - Casimir like forces CHU, after the Dutch physicist 
Hendrik Casimir who in 1948 predicted an attractive force be¬ 
tween two parallel perfectly conducting metal plates (T) sepa¬ 
rated by a finite gap L in vacuum at zero temperature. In order 
the force to be long-ranged (i.e., to decay in a power-law and 
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not in an exponential with the distance way) the interactions in 
the system have to be mediated by massless excitations - pho¬ 
tons, Goldstone bosons, acoustic phonons, etc. Considered 
in this general form the Casimir effect is a subject of investi¬ 
gations in condensed matter physics, quantum electrodynam¬ 
ics, quantum chromodynamics and cosmology. The results 
are summarized in an impressive number of reviews IMS). 

When the fluctuating field is the electromagnetic one the 
effect is known as the quantum electro-dynamical (QED) 
Casimir effect. There the Casimir force is caused by zero- 
point and thermal fluctuations of the electromagnetic field. In 
first approximation, it depends only on the velocity of light 
c, Planck’s constant Tt, the temperature T, and the separation 
distance between the bodies L, i.e., this force to a great extend 
is universal. A more advanced theory, the so-called Lifshitz 
theory, reveals the dependence of the Casimir force on the ma¬ 
terial properties of the bodies l35lj37l l and geometry of their 
boundary surfaces l38l[39j . 

Thirty years after Casimir’s prediction, M. E. Fisher and P. 
G. de Gennes suggested that the fluctuating medium, confined 
between the bodies can be a fluid, the fluctuating field be¬ 
ing the field of its order parameter, in which the interactions 
in the system are mediated not by photons but by different 
type of massless excitations like critical fluctuations or Gold- 
stone bosons (spin waves). The corresponding Casimir effect 
is known as the thermodynamic Casimir effect 0- When the 
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confined fluid approaches it’s critical point, the correspond¬ 
ing fluctuations are the critical fluctuations of the order pa¬ 
rameter and then the effect is usually called critical Casimir 
effect. In first approximation, the thermodynamic Casimir ef¬ 
fect depends only on the gross features of the system - its di¬ 
mensionality d and the symmetry of the ordered state n (both 
defining the so-called bulk universality class of the system) 
and on the boundary conditions (determined by the surface 
universality classes). Therefore, to a great extend the thermo¬ 
dynamic Casimir force is also universal. So far the critical 
Casimir effect has enjoyed two general reviews ll28i j29j and 
some concerning specific aspects of it 1301434) . 

Currently the Casimir effect is an object of intensive stud¬ 
ies both in its original formulation due to Casimir as well as 
in its thermodynamic manifestation. In the present article we 
will report theoretical results dealing with the critical Casimir 
effect. Let us note, that the critical Casimir effect has been 
already directly observed, utilizing light scattering measure¬ 
ments, in the interaction of a colloid spherical particle with 
a plate ll42l both of which are immersed in a binary liquid 
mixture. The effect has been also studied in 4 He 8431 . 1441 . 
as well as in 3 He- 4 He mixtures m in the context of forces 
that determine the properties of a film of a substance in the 
vicinity of its bulk critical point. In Ref. 8461 and Ref. 147) 
one has performed measurements of the Casimir force in thin 
wetting films of binary liquid mixture. One the theoretical 
side, the effect has been studied via exact calculations in the 
two-dimensional Ising model EEH52, the three dimensional 
spherical model 15814661 . via conformal-theoretical methods 
EZHEl, within mean-field type calculations on Ising type 
Gma and XY models (80), through renormalization-group 
studies via £-expansion MM and via fixed dimension d 
techniques 8871489) of 0(n) models, as well as via Monte- 
Carlo calculations (791 190lfl03l. The fluctuation of impor¬ 
tance in all of the above mentioned models are of thermal 
origin since all these models possess non-zero critical tem¬ 
perature. In some systems, however, certain quantum param¬ 
eters govern the fluctuations near their critical point which is 
usually close to or at the zero temperature 1104141071 . In this 
particular case one speaks of a quantum critical Casimir effect 
829l lK)8l[T09l . 

The rapid progress in nanotechnology has resulted in the 
growth of interest in fluctuation-induced phenomena, which 
play a dominant role between neutral non-magnetic objects 
at short separation distances (below a micrometer). The van 
der Waals and Casimir forces, both known under the generic 
name dispersion forces, play a key role in Micro- and Nano- 
Electromechanical Systems (MEMS/NEMS) 111044112l oper¬ 
ating at such distances. Indeed, upon scaling down devices, 
the dispersion forces can induce some usually undesirable 
non-linear behaviors in such systems cm Irreversible phe¬ 
nomena appear such as stiction (i.e., irreversible adhesion) 
or pull-in due to mechanical instabilities 11 1411115 1. There¬ 
fore the ability to modify the Casimir interaction can strongly 
influence the development of MEMS/NEMS. Several theo¬ 
rems seriously limit, however, the possible search of repul¬ 
sive Casimir forces (iMm). Currently, apart from some 
suggestions for achieving Casimir repulsion in systems out of 
equilibrium ll2Tl|ll9l4l27ll . the only experimentally well ver¬ 
ified way to obtain repulsive Casimir force is to have interac¬ 
tion between two different materials characterized by dielec¬ 


tric permittivities £i and £2 such that EH37) 

£l<£ M <£2 (1.1) 

along the imaginary frequency axis, with £ w being the dielec¬ 
tric permittivity of the medium in between them. In Refs. 
Il28lll33l Casimir repulsion was indeed observed experimen¬ 
tally for the sphere-plate geometry. 

In the current article we study the interplay between the 
critical Casimir force and the van der Waals one in a system 
composed out of two flat parallel slabs both immersed in a 
critical fluid. Let us note that both the critical Casimir and 
van der Waals forces are fluctuation induced ones but due to 
the fluctuations of different entities. For terminological clar¬ 
ity, let us also remind that in colloid sciences fluid mediated 
interactions between two surfaces or large particles are usu¬ 
ally referred to as solvation forces |48l [134 , 1351 . Thus, we 
study here a particular case of such a force when the fluid 
is near its critical point. In our system we suppose that the 
slabs are coated by thin layers of some substances, confining 
either a non-polar one-component fluid or a non-polar binary 
liquid mixture. We suppose that the liquid phase of the one- 
component fluid or one of the components of the binary liquid 
mixture are strongly adsorbed by both coating layers, i.e., they 
ensure the so-called (+,+) boundary conditions. The slabs, 
on the other hand, influence the fluid by long-range compet¬ 
ing dispersion potentials. In the case of a simple fluid these 
potentials increase the adsorption of the fluid near one of the 
surfaces, leading to preference there of its liquid phase, and 
decrease it near the other one. In the case of a binary fluid 
mixture the substrates prefer one of the components near the 
top and the other one near the bottom of the system. We 
will demonstrate that this experimentally realizable compe¬ 
tition between the effects due to the coatings and the slabs can 
result in interesting effects like sign change of the Casimir 
force, acting between the surfaces confining the fluid when 
one changes t, ll or L. The last facts can potentially be used 
in designing nanodevices and for governing the behavior of 
objects at small, below micrometer, distances. 

If a fluid system possesses a surface it breaks the spatial 
symmetry of the bulk system. The quantitative effects of the 
presence of a surface on the thermodynamic behavior of the 
system depends on the penetration depth of this symmetry 
breaking effect into the volume. There are two phenomena 
which increase the surface effects: long-range interactions 
and long-range correlations. They can act separately, or si¬ 
multaneously, which leads to an interesting interplay of the 
effects due to any of them ll76l 113644139) . The penetration 
depth due to the correlations is set by the correlation length 
(j of the order parameter of the system; becomes large, 
and theoretically diverges, in the vicinity of the bulk critical 
point (T c ,p c ): Z{T^T+,n= p c ) - t = (T- T c )/T c , 

and t, (T = T Cl ji ->• jx c ) ~ £o„u \Ap/{k B T c ) |“ V / A , Ap = p - p c , 
where v and A are the usual critical exponents. If the system 
is made finite, e.g., by the introduction of a second wall the 
behavior of the fluid is further enriched. When E, becomes 
comparable to the characteristic system size, say L, the size 
dependence of thermodynamic functions enters through the 
ratio L/^, i.e., takes a scaling form given by the finite-size 
scaling theory lf29l 40] 4T] jl40|4l42l that incorporates, inter 
alia, shift of the critical point of the system 1143141471 . Be- 
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low T c if the confining walls of the film geometry consist of 
the same material, one encounters the phenomenon of capil¬ 
lary condensation f76l 134 . 14711149 1 where the confinement 
of the fluid causes, e.g., the liquid vapor coexistence line to 
shift away from the coexistence line of the bulk fluid into the 
one-phase regime. We will demonstrate in the current article 
that in the envisaged realization of our system depending on 
the material properties of the slabs the phase diagram of the 
finite system might essentially differ from that one of the well 
studied case of a system with short-ranged type interactions 
and strongly adsorbing surfaces. 


The article is arranged as follows. In Sec. [0]we recall and 
comment on the finite-size behaviour of systems with disper¬ 
sion forces extending the known facts to the expected behavior 
of the Casimir and the net forces when they act between walls 
being semi-infinite slabs coated by some thin substances. By 
doing so, we especially pay attention to the conditions under 
which the effects stemming from these interactions are rele¬ 
vant. Section |III] presents the corresponding lattice gas mod¬ 
els suitable for the investigation of fluid media with account 
of the long-ranged van der Waals interactions. Here we iden¬ 
tify the main coupling parameters characterizing the interac¬ 
tions in the systems and in Sec. IV the equation for the equi¬ 
librium profile of the finite-size order parameter is obtained, 
which we later use to calculate the forces of interest. Section 
[VJpresents the numerical results for the behaviour of the inves¬ 
tigated forces followed by Sec. VI where the phase behaviour 
of the considered type fluid system is briefly discussed. The 
experimental feasibility of the predicted effects is discussed 
in Sec. VII Here we also comment on the possible applica¬ 
tion of our findings in design of nanodevices and for govern¬ 
ing the behavior of objects, say colloidal particles, at small 
distances. The article ends with a summary and discussion 
section - Sec. VIII Important technical details concerning 
the Hamaker term for a van der Waals system of two differ¬ 
ent substances separated by a fluid medium are presented in 
Appendix [A] 


II. THE THERMODYNAMIC CASIMIR FORCE IN A 
NON-POLAR FLUID FILM SYSTEMS WITH DISPERSION 
FORCES 


Let us consider some fluid medium confined between two 
parallel slabs of some materials Si and S 2 . Any of the slabs 
is coated by thin solid films of some other substances Li and 
L 2 , respectively (see Fig. [TJ. Let the slabs are situated at some 
distance L from each other. We suppose that the thicknesses 
of the coating films are negligible. In the remainder of the text 
we are going to designate each slab and the thin solid film that 
coats it as a "wall", and refer to any of its two components 
separately only when this is necessary. 

If the fluid medium is in contact with a particle reser¬ 
voir with a chemical potential jj., the grand canonical poten- 
tial Q} ex (L\T,n) of this medium in excess to its bulk value 
g/LcObuikiT, fi) depends on L and, thus, one can define the ef- 
fective force F^ ol (L\T,fj.) per cross sectional area .t/ and k/fl\ 
due to the fluctuations of the medium and dispersion interac- 






of liquid B molecules 


FIG. 1. Schematic depiction of a finite-size fluid system, consist¬ 
ing of two parallel slabs of some substances §1 and §2> coated by 
thin layers of some other substances Li and L 2 , respectively, confin¬ 
ing some fluid medium M - a non-polar one-component fluid [as an 
example we depicted the carbon dioxide molecules (CO 2 )] (a) or a 
binary mixture (b) composed out of the molecules of the non-polar 
liquids A and B [the depicted examples include tetrachloroethene 
(C 2 CI 4 ) - as the substance A, and benzene (CgHg) - as the sub¬ 
stance B\. The confined fluid medium is considered embedded on a 
lattice in which (a), some nodes are occupied by a particle and others 
are not - thus depicting the "liquid" and "gas" states respectively at 
some values of fluid temperature T and chemical potential fl, or (b) 
some of the nodes are occupied by a molecule from the substance 
A (the "liquid" state) and the rest are occupied by the molecules be¬ 
longing to the species B ("gas" state). The confining walls impose 
on the fluid medium boundary conditions of strong adsorption on the 
coating layers, i.e., the nearest to the coating substances layers are 
entirely occupied by the particles of the one-component fluid or if 
the medium is a binary liquid mixture - by the particles of one of its 
components (in the presented figure we choose the molecules of the 
specie A). 




tions in it as 

pF^(L\T,n) =f${L\T,n) = (2 1} 

where the superscript T designates the boundary conditions 
which the confining walls impose on the fluid medium 
(see above), (o£\l\T,h) = <o^(L\T,n) ~ = 

Q.il\L\T,fl)/&/ is the excess grand canonical potential per 
unit area g/, (L\T.ji) = srf fflW (L\T, fi) is the total grand 

canonical potential, fflbuik(T,/i) ' s the density of the bulk 
grand canonical potential, and p = 1 /(ksT) 111501 . Let us 
stress that, as pointed out in Ref. It76ll . one should keep in 
mind that the force f tot (L\T, fl), see Eq. m depends on how 
one defines the thickness of the film. This implies that a quan¬ 
titative comparison between experimental data and theory is 
only possible if the data are accompanied by a precise defini- 
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tion of what L is. 

Away from the critical temperature of the system it is cus¬ 
tomary to write the force acting between the plates of the fluid 
system in the form 

f tot {L\T,n) ^ (<7 - l)l5H A (T,n)L- a ^ d , ( 2 . 2 ) 

where one normally considers the case d = a and omits the 
apparent dependence on the so-called retardation length [143 , 

1151 1 ^ re t- Here H A is the Hamaker term, whose dependence 
from the temperature and chemical potential is given by the 
so-called Hamaker constant If54l[l52l 

A Uam {T,n) = -\2nH A (T ,n). (2.3) 

The Hamaker constant, as it is clear from above, is a con¬ 
stant only in the sense that it is L-independent. It depends on 
the temperature, chemical potential and on the material prop¬ 
erties of the fluid medium and the walls. The factor 1 2it in 
Eq. ( |2.3| ) is introduced there due to historical reason, accord¬ 
ing to which the interaction energy between two substrates at 
a finite separation L in the case of standard van der Waals in¬ 
teraction (i.e., d = (7 = 3), away from any phase transition 
region, isEHna 

co ex (L\T,ii) = ~-^-A mm (T,ii)L- 2 . (2.4) 

127T 

The Hamaker term takes into account the leading L-dependent 
parts of the i) direct interaction between the slabs A Sl)S2 , ii) 
between each slab and the fluid medium - A Sli /, and A S2j /, as 
well as the Hi) interactions between the portion of constituents 
of the fluid medium situated within the cavity bounded by the 
substrates - A/, i.e., 

H A (T,h) =A l (T,n)+A SltS2 (T) 

+A Slt i(T,n)+A S2t i(T,n). (2.5) 

Note that in Eqs. d 22) - £5) both the slabs and the fluid 
medium are characterized by their bulk properties at the given 
temperature and chemical potentials. 

Near the critical temperature T c of the bulk system Eq. \2.2\ 
is no longer valid since the critical fluctuations of the order 
parameter lead to new contribution to the total force called 
thermodynamic (critical) Casimir force (see below). For such 
a system, following Ref. EE near the bulk critical point 
Eq.dTT} can be written in the form 

/tot(L|r,/x)~L -Xcrit [x t ,Xn,Xi,{x Sn i= 1 , 2 } 

+{a-\)^H A {T^)L-°^- d . ( 2 . 6 ) 

In Eq. \2.6) X cr -,i is dimensionless, universal scaling func¬ 
tion, x t = t (L/^ 0 h ) 1// ' and = f}A)X (L/^o^) 4 ^ are the 
temperature and field relevant scaling variables, respectively, 
while xi = A(L/£q“) m ‘, x Si = s-, (L/£,q ) Ws , i = 1,2 and 
x g = gai (L/£ 0 + ) m are the irrelevant in renormalization group 
sense scaling variables associated with the interactions in the 
system. The freedom of choosing the precise definition of 
what L is in systems with boundaries leads to the formal ne¬ 
cessity to write L as L + Lq, with Lq being a microscopic 
length. That will lead to further scaling corrections propor¬ 
tional to L _1 . Since in the current article we will keep in all 


calculated quantities only their leading L dependence, we re¬ 
frain from further refinement of the scaling Ansatz ( |2.6| >. By 
comparing Eq. \2.2\ and Eq. ( | 2.6 [ > one immediately concludes 
that A cr j t tends to zero away from the critical point, i.e., when 
at least one of the relevant scaling parameters \x t \ and \x^\ 
becomes large, i.e., when \x t \ 3> 1 and/or |jq, | 1. 

As it is well known, the critical behavior of simple fluids 
and of binary liquid mixtures is described within the Ising uni¬ 
versality class which determines the values of the critical ex¬ 
ponents A = j3 + y and v ones the dimensionality of the system 
d is fixed. When d = 3 this universality class is characterized 
by critical exponents ll 1531 

v = 0.630, J3 = 0.327, y = 1.237,0 = GJv = 0.524. (2.7) 

In order to better reflect the actual properties of the non-polar 
fluids, instead of considering nearest-neighbour interactions 
we assume long-ranged pair ones between the fluid particles, 
decaying asymptotically ~ J l r~ d a for distances r between 
each other, and substrate potentials ~ J' , j z2 a ■ i = 1,2 acting 
on the fluid particles at a distance z from each of the two slabs. 
We recall that when a >2 systems governed by such long- 
range interactions, usually termed subleading long-ranged in¬ 
teractions nmum also belong to the Ising universality 
class characterized by short-ranged forces lll56l . The last im¬ 
plies, among the others, that the critical exponents, e.g., do 
not depend on a for such type of interactions. An important 
representative of such type of interactions are the non-retarded 
dispersion interactions with d = <7 = 3, which are one of the 
three types of van der Walls interactions. By varying the ratio 
between the strengths of the long-ranged - J 1 and the short- 
ranged - J l SI contributions one can quantitatively probe the 
importance of the long-ranged parts of the interactions within 
the fluid medium and study potential experiments in colloidal 
systems which allow for a dedicated tailoring of the form of 
the effective interactions between colloidal particles. 

In Eq. ( |2.6| l, <37 is the standard correction-to-scaling expo¬ 
nent for short-range systems, while 07/ = a — (2 — T )) and 
GJ S = a — (d + 2 — rj)/2 are the correction-to-scaling expo¬ 
nents due to the long-range parts of the interaction potentials 
between the constituents of the fluid medium and those of the 
confining walls. Further L—dependent contributions to the 
total forces / tot such as next-to-leading order contributions 
to the Hamaker terms or higher order corrections to scaling 
are neglected because they are smaller than those captured in 
Eq. ( |2.6| ). The exponent tj, which appears in the expressions 
for GJ/ and E7 V , is the standard one characterizing the decay of 
the bulk two-point correlation function at the critical temper¬ 
ature, g m is the (dimensionless) scaling field associated with 
the Wegner-type corrections, while A and s,, i = 1,2 are di¬ 
mensionless non-universal coupling constants: A is propor¬ 
tional to the strength J 1 of the long-range part of the inter¬ 
action potential between the particles of the fluid, whereas 
Si, i = 1,2 are proportional to the contrast between the po¬ 
tentials of the bounding slabs and those in the fluid medium 
(see below). For systems belonging to three-dimensional Ising 
universality class with "genuine" non-retarded van der Waals 
interaction one has d = cr = 3, and t] = 0.03627(10) lll53l . 
This leads to GT, ~ 1.03, Q s ~ 0.52, and G7 = 0.832(6) fl53l . 
Within the mean-field theory with d = a = 4 and r) = 0 one 
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has, instead, 01/ mf = 2 and CJ^mf = 1 ■ One then has 

/ tot (L|r, J u)L 4 ~X crit (-) + 3/3// A (7’, J u), (2.8) 

where all critical exponents take their mean-field values /3 = 
V = 1/2, 7=1, A = 3/2. 

The peculiarities of the scaling theory for systems with 
dispersion forces are described in Refs. IhTl [76l 113611137 , 
Il54l 1155 1. One obtains that despite these systems do be¬ 
long to the Ising universality class with short-ranged forces, 
the finite-size quantities decay algebraically with L towards 
their bulk values, and not in an exponential in L way, when 
4 ln(^/ao) L even for ao, where ao is the charac¬ 

teristic distance between the molecules of the fluid system. In 
this regime the dominant finite-size contributions to the free 
energy and to the force between the walls bounding the sys¬ 
tem stem from the long-ranged algebraically decaying parts of 
the interaction potentials. One can formulate a criterion clar¬ 
ifying when the long-ranged tails of the interactions can not 
be disregarded even in the critical region of the finite system. 
Neglecting the thickness of the coating layers, for the system 
under consideration the corresponding criterion states that the 
long range tails of the interactions can be disregarded only 
when 1611 [761 fl37l 

2 ff (| Sl | + | S2 |)[L/^] A/V “ CT <l. (2.9) 

Using the scaling relations, the above can be formulated as a 
constraint on the thickness of the system under study. One 
obtains that when L L cr it, where 


rections to the leading L dependance. The reader can refer for 
a detailed comment on that matter to Ref. Il76l (see also Ref. 
ED for the properties of X^). The validity of the proposed 
criterion as well as the statements made beneath it are well 
illustrated on Fig. [2] 

The contribution of the dispersion forces to the total effec¬ 
tive force / tot can be distinguished from that of the critical 
Casimir force by their temperature dependence, because the 
leading temperature dependence of the former does not ex¬ 
hibit a singularity. Thus, one has 

f tat (L\T,n)=f^(L\T,n)+f^ s) (L\T,n), (2.13) 

where 

f vm = f^ s \L\T,^), (2.14) 

and 

f5F A ^(L\T,n)=f^ s \L\T,^). (2.15) 

One expects that near the bulk critical point 

PF AtQa (L\T,n) =L- d X Cas (x t , Xtl ,---), (2.16) 


where A/as is a scaling function that for large enough L (see 
below) approaches the scaling function of the short-ranged 
system A/ r as (x t ,x^). From Eqs. (2.6 1 , (2.13 i - (2.16 1 it fol¬ 
lows that the scaling function of the critical Casimir force A/as 
is proportional to the sum of X cr it and the singular part of the 
Hamaker term. 


L c rit = W[2 tT (N + M)] v//3 I (2.10) 


the effects of the long-ranged tails of the interaction can be 
neglected within the critical region of the finite system. For 
d = c7 = 3 and with /I and v from Eq. for the three- 
dimensional Ising model one obtains 

L crit ~ 55£o,//(M + M) 1 ' 923 . (2.11) 


We conclude that for moderate values of L, i.e., when L < 
L cr it the behavior of the Casimir force can strongly depend on 
the details of the fluid-fluid [the amplitude ^o,/i depends on 
the details of the fluid-fluid interaction (see Eqs. (4.15) and 
(4.17) in Ref. |[76l and the text therein)] and substrate-fluid 
interactions which, as it turns out, can influence even the sign 
of the force. For large L, i.e., for L 3> L cr it the behavior of the 
force shall approach the one of the short-ranged system. The 
last can be easily seen from Eq. ( |2.6| ) by simply expanding 
then the scaling function A^ c rit 

/tot(T|L,/x) — L |^cri t ( x «J X ju) +X/Xcrit(^/,X jU ) 

+ \ [ x siKk( X t^n)+ X S2 X ^t( x t^n)] j 

+(C7- 1)PH A (T, (2.12) 


Here A// originates from the short-range interactions [see 
Eqs. (4.10 1 ], It is well known that X^ t < 0, i.e., the force 
is attractive, under (+, +) boundary conditions. A/ r it provides 
the leading behaviour of the force near the bulk critical point 
(x, = 0, x M = 0). There X' m . X^' a and A// represent only cor¬ 


We will often compare the behavior of the system with sub¬ 
leading long-ranged interactions present with this one of a sys¬ 
tem with purely short-ranged interactions which will serve as 
a reference system. In such a purely short-range system one 
has H a = 0. Then, at the bulk critical point (T = T c , /i = fi c ) 
the leading term of the thermodynamic Casimir force between 
the slabs bounding the fluid has the form 

= (d- l)A«(J)^J, (2.17) 

where X*^ t (0) = (d — l)A^(of). Here A ^\d) = &{\) is an 
universal dimensionless quantity, called Casimir amplitude, 
which depends on the bulk and surface universality classes 
(specified by the boundary conditions t). Since the Casimir 
force is proportional to kgT c the interaction between the walls 
can become rather strong in a system with high critical tem¬ 
perature such as, e.g., in classical binary liquid mixtures. Note 
that the sign of the force depends on the sign of the Casimir 
amplitude A ir> (d) which, on its turn, depends on the boundary 
conditions t. According to the usual convention negative sign 
corresponds to attraction, while positive sign means repulsion 
of the surfaces bounding the system. 

The experimental and theoretical evidences accumulated 
till nowadays support the statement that the Casimir force is 
attractive when the boundary conditions on both plates are the 
same, or similar, and is repulsive when they essentially differ 
from each other, e.g., when in the case of a one-component 
fluid one of the surfaces adsorbs the liquid phase of the fluid 
while the other prefers the vapor phase. It is instructive to go 
back to the explicit physical units in Eq. \2.\1\ for the physi- 
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* ( =t(L /£ 0 + ) 2 x„=p e Ap(L/^o ,„) 3 




x^PcAniL/go^) 3 


FIG. 2. Behaviour of the scaling function X CI it{xt,Xn) in ad = 3 dimensional confined fluid system. On (a) and (b) the parameters reflecting 
the substrate-fluid interactions have values V| c =1.0 and si tC = —1-0, while on (c) and (d) we took s\ c = s 2 c = —1-0. In all four cases the 
parameter reflecting the role of the long-ranged fluid-fluid interaction is A = 2.0. The considered separations between the confining walls are 

L = 20 (“■“), L = 90 (-) and L = 800 (“♦“). In (a) and (c) the functional dependence on the temperature scaling variable x, at Afl = 0 

is visualized while in (b) and (d) the dependence on the field one Xn at T = T C is presented. As it can be seen for L L cr ; t — 90, the van 
der Waals interactions between the walls and the fluid medium become irrelevant and the scaling function for a fluid system with L = 800 
(-♦-) is very similar in shape to a short-range one (“•”). Due to the positive substrate-fluid parameter in cases (a) and (b) this effect is 
more pronounced than in (c) and (d). For smaller separations L = 20 (“■”) the effect of the long-range van der Waals interactions becomes 
relevant everywhere, including the critical region of the system where the scaling function changes it’s sign twice [see (b) and (c)] or even 
three times [see (d)]. In case (a) we observe that while for L = 20 (-■-) as a function of x, at Xn = 0 the corresponding force is everywhere 
repulsive, it becomes under enlarging L everywhere attractive, say, for L = 800 ("♦“), as in the case of a system with completely short-ranged 
interactions (“•-). When fl = jl c the curve for the short-ranged interactions depicts the numerical evaluation of the exact analytical result 
given by Eq. ([4. 1 0), while for T = T C this curve is evaluated within the presented mean-field theory. 


cally most relevant case of d = 3. One has 


under which our results have been obtained. 


i&a)*..ix.o- AW <"= 3 > 4 N 


(L/iim) 


(2.18) 


where T mom = 20 °C (293.15 K). Sine as discussed above, for 
most systems and boundary conditions A ^(d) = when 
T c r — T room the thermodynamic Casimir force, for some space 
separation L, will be of the same order of magnitude as the 
quantum one 


cQED 
F A.C as 


(L) 


! 3x10-3 1 N 

240 (L/ fim) 4 m 2 


(2.19) 


and they both shall be significant and consequently measur¬ 
able at or below the micrometer length scale. 

We turn now to description of the model and the procedure 


III. THE MODEL 


We are going to utilize the same type of model already used 
in Refs. Il76l 1 136 1 but amended to take into account the spe¬ 
cific features of the system considered in the current article. 
Among them is the role of the two competing substrate poten¬ 
tials. Here, in order to introduce the notations needed further, 
we briefly recall the basic expressions of that model paying a 
bit more attention only to difference of the current model with 
that one studied in Refs. (761 11361 . 

We consider a lattice-gas model of a fluid confined between 
two planar walls, separated at a distance L from each other. 
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with grand canonical potential £1 [p(r)] given by 

Q.[p(r)]=k B T £ {p(r)ln[p(r)] 

reM 

+ [l-p(r)]ln[l-p(r)]} + i £ p(r)tv'(r-r')p(r') 

Z r,r'GM 


+ £ M s >l^)(z)-plp(r), 


(3.1) 


where M is a simple cubic lattice in the region occupied by 
the fluid medium - oo^ 1 x [0,L] and W 1 * (z) is an external 

potential that reflects the interactions between the confining 
walls and the constituents of the fluid, given by 

yN^) (z) = -p Sl J^5(z)-p S2 J s i’ l 8(L-z) 

+v Sl (z+l) a + v S2 (L + 1 — z) a , (3.2) 

where v Sj = —G(d, o)p Sj J Si ’ 1 , i = 1,2, with 

F () 

G(d,o) =AK {d - 1 ^ 2 — 1 ? 4 . (3.3) 

V 7 aT{^) 


This type of functional can be viewed as a modification of 
the model utilized by Fisher and Nakanishi 01431 1 1571 in their 
mean-field investigation of systems governed by short-range 
forces. 

In Eq. m the terms in curly brackets multiplied by kgT 
correspond to the entropy contributions to the total energy, 
w l (r — r') = —4/ , (r —r') is the non-local coupling (interaction 
potential) between the constituents of the confined medium 
and p is the chemical potential. Here and in the remainder 
of this paper, all length scales are taken in units of the lattice 
constant ao (for concrete values see Table. [I}, so that the par¬ 
ticle number density p (r) becomes simply a number density 
which varies in the range [0,1], 

The variation of Eq. with respect to p (r) leads to the 
equation of state for the equilibrium density p*(r) 


2p*( r ) 


= tanh 


f E w'(r —r')p*(r') 

Z r'eM 


4 


'n-vMM(zj 


(3.4) 


The advantage of this type of equation is that it lends itself to 
numerical solution by iterative procedures. For a given geom¬ 
etry and external walls-fluid potential y( 11 l / l i2 )(^) its solution 
determines the equilibrium order parameter profile p*(r) in 
the system. Inserting this profile into Eq. ( |3. 1 1 >, renders the 
grand canonical potential of the considered system. 

Denoting, as in Refs. f76l [136 1 <j>*( r) = 2p*(r) — 1 and 
Ap = p — Pc, where p r = ^ £ r < vv ; (r — r'), the equation of state 
Eq. (|3.4|> can be rewritten in the standard form 


0*(r) = tanh< (3 £/( r-r>V) 

[ r'eM 

+ |[Ap-Ay(z)]|. (3.5) 

The bulk properties of the model are well known (see, e.g. 
1 158] [159 1 and the references therein). We recall that the or¬ 


der parameter p* of the system has a critical value p* = 0 
which corresponds to p c = 1/2 so that p* = 2(p* — p c ). 
The bulk critical point of the model is given by {/3 = j3 c = 
[Lr ( r )] -1 > P = Pc = —2 £ r f (r)} with the sum running over 
the whole lattice. Within the mean-field approximation the 
critical exponents for the order parameter and the compress¬ 
ibility are fi = 1/2 and 7=1, respectively. The effective sur¬ 
face potential Ay (z) in Eq. (|3.5|> is given by 


Ay(z) 


<5v S| 5 v S2 

(z+l) CT + (L+ 1 —z) a 1 


(3.6) 


where 1 < z < L — 1, and contributions of the order of z a 1 
z^ CT_2 , etc., have been neglected. 


<5v. Si - -G(d, a) {p s f- 1 - p c J ') ,i = 1,2 (3.7) 


are (T- and p-independent) constants, 


/(r)=4{5(|r|) + 5(|r|-l)} + 


M H-i) 

l + |r|rf+u ’ 


(3.8) 


is a proper lattice version of — w ? (r)/4 as the interaction en¬ 
ergy between the fluid particles, and 


/^(r)=/2-'5(|r|-l) + 


J Si ' l e (|r| — 1) 


• \d+a 


, * = 1,2 (3.9) 


is the one between a fluid particle and a substrate particle, 8 (x) 
is the discrete delta function and 6(x) is the Heaviside step 
function with the convention 0 (0) = 0; in Eq. ( | 3.9 [ > p Sj , (=1,2 
are the number densities of the coated slabs in units of d (for 
concrete values see Table. 0- Note that the effective poten¬ 
tials (Siy, i = 1,2 result from the difference between the rel¬ 
ative strength of the substrate-fluid interactions for substrates 
with density p Si , (=1,2 and that of the fluid-fluid interactions 
for a fluid with density p c . In Eq. ( |3 .6 [ > the restriction z > 1 
holds because we consider the layers closest to the substrate 
to be completely occupied by the liquid phase of the fluid 
(which implies that we consider the strong adsorption limit), 
i.e., p(0) = p{L) = 1, which is achieved by taking the limits 
J%' 1 —► °°, i = 1,2; thus, the actual values of Ay(0) = A V(L) 
will play no role. In order to preserve the monotonic behavior 
of vv ; (r) as a function of the distance r between the particles, 
in Eq. ( |3.8| ) we have to require that j\ x > J l /(\ + 2 d+a ), i.e., 
A < 1 +2 d+tT , where 


A =i'/4- (3-10) 

From Eqs. ( |3.5| > - ( | 3.9 [ i one can identify the dimensionless 
coupling constants y, (=1,2 appearing in Eq. p.6| l 

Si = ~p8v sn (=1,2. (3.11) 

Here ,v, > 0, i.e., p Sj J Si ’ 1 > p c J 1 corresponds to walls "prefer¬ 
ring" the liquid phase of the fluid, while ,v, < 0, or p Si J s ’’ 1 < 
PcJ 1 mirrors the one with affinity to its gas phase. The 
marginal case s,- = 0 which corresponds to p Si J suI = p c J l will 
be commented further in the text. 


In terms of p the functional defined by Eq. m takes the 
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form 


«[ 0 (r )]=k B T£ 

reM 


n+^in 

[l + 0(r)l 

l 2 

2 


+ ‘-» (r > i„ 
2 


l-0(r) 

2 


--£[A M -AV( z )]0(r) 


' reM 


E ^(r-r')0(r)0(r')+a e g, 


(3.12) 


where 


^ re g — 9 


reM 


AH-A V(z)~ E ^ Z (r-r') 

r'eM 


(3.13) 


do not depend on 0 (r) and, therefore, is a regular background 
term carrying an L —dependence and thus showing up in the 
corresponding force. 

The only quantity in Eq. ( |2.6| l, which still has to be iden¬ 
tified for our model, is the value of the coupling constant A. 
According to Ref. m 


A = v a /v 2 , (3.14) 

where v a and V 2 are coefficients in the Fourier trans¬ 
form J\ k) = £ r exp(/k • r)J l {r) of the interaction J l ( r) [see 
Eq. ( |3.8| >] (for more details see Ref. Ii76l ). It turns out that A 
depends on A. 

In accordance with Eq. for the finite-size behavior 

of the excess grand canonical potential per unit area srf of a 
liquid film in the case when both confining surfaces strongly 
adsorb the liquid phase, one expects 

ft) ex (L|r, M )^yr + )£ s + ^An(*) 

+H A (T,n)L-^^ d , (3.15) 

where yj ls and $ ls are the non-singular parts of the surface ten¬ 
sions at the surfaces of the confining walls, while the singular 
parts are incorporated in the scaling function Xq(*), which 
arguments are as those of the function A cr j t in Eq. ( |2.6[ >. 

In Eq. ( |3.15| > H A is the Hamaker term which explicit form 
is derived in Appendix [A] The result is 

WS) = ( P “ - T^ n ) ( p “ - £»)' (316) 

Note that this result is in full agreement with the 
Dzyaloshinskii-Lifshitz-Pitaevskii (DLP) theory l35l l37l . It 
provides, however, an easy possibility to study the T and p 
dependence of H A by studying the corresponding dependen¬ 
cies of pt,, say, near the critical point of the fluid system. In 
accord with the DLP theory it teaches us that the sign of the 
Hamaker term depends on the contrast of material properties 
of the two bounding substrates with respect to the substance 
which is in-between them. Coating the substrate surfaces of 
the system with some additional material does not change the 
leading-order L dependence of the interactions between the 
substrates and, therefore, does not change the above property. 
Let us stress that when si = S 2 = s, i.e., when the bounding 
substances are made from same material, the Hamaker term 
will be negative, independent on the properties of the fluid 


in-between them, i.e., 


H a = —C(d,o)J Sl,Si 



i. ) 


2 


< 0 , 


(3.17) 


which corresponds to attraction between the confining walls. 
Similar is the situation when pb —> 0, i.e., when the fluid sep¬ 
arating the substrates is replaced by vacuum. Then, again, 
H a < 0 irrespective of the material properties of the bounding 
substances, since linip^ a)H,\ = 4 S| V2 < 0. Further details on 
the behavior of Ha are given in Appendix |A| 

In the next section IV we are going to present results for 
(O ex (L\T,p) based on the model described in the current sec¬ 
tion. 


IV. FINITE-SIZE BEHAVIOUR OF THE MODEL IN A 
FILM GEOMETRY 


Let us start by rewriting the equations presented in Sec. Ill 


in the form suitable for studying of a system in a film ge¬ 
ometry. Because of the translational symmetry of the system 
along the bounding surfaces, the quantities of interest depend 
only on the spacial coordinate z. along which the system is 
finite. Thus, one can write: 0(r) = <j> (r||,z) = 0(z), where 
r = {r||,z}, i.e., the local order parameter profile is given by 
{0(z), 0 < z < L}, with 0(0) = 0(L) = 1. Hence the equation 
Eq. (|3.5|> for the equilibrium profile becomes 


arctanh[0*(z)] = ^ [Aju — AV(z)] 

+K { aj a (A) r (z) + (A) [0*(z +1)+ 0*(z-1)] 

E ft/.aflz-z'IMIz-z'l - WV)} , (4-1) 

z'=0 J 


where Aju = p — p c , AV (z) is defined in Eq. ( |3.6| >, K = pj' sr , 
a d , a (A) = {2d - 1) + A (c d , a - d) , aj a (A) = 1.0 + A {cj a - 
0.5) with cj a = g di a(l) +g“ ff (±l)- The functions c d ^ a , 
gd,o(\z-?!\) and g n da {\z~z!\) are determined in Eqs. (CIO), 
(Cl 1) and (C12) of Ref. Il76l . respectively. 

When L —> °° the equilibrium finite-size order parameter 
0*(z) tends to its bulk value 0/, = 2 (pb — p c ) € [0,1]. The 
corresponding equation for &/,, performing the limit L —> °° in 
Eq. (|4. 1 [>. reads 


(pb = tanh 



+K 


a d ,a (A) + 2 a d " a (A) + 2 A E 8d,o{z ) 

J>2 



(4.2) 


wherefrom one immediately identifies the coordinates of bulk 
critical point 


K c { A) 


&d,o (A) + 2 a™' c (A) + 2A E gd,a{z) 

z>2 


(4.3) 


and Ap = 0. Note that the position of the bulk critical point 
depends on A, i.e., on the presence and the strength of the 
long-ranged tails in the fluid-fluid interactions. 
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For a fluid confined to a film geometry the natural quantity 
to consider is the excess grand canonical potential per unit 
area stf\ co ex = lirn K /_ >00 — LObuik- Using the result 

of Ref. [132] [see Eq. (3.14) there], as well as the identity 
arctanh 0 = ( 1 / 2 )[ln( 1 + 0 ) — ln(l — (j))], one can write j 8 ® ex 
in the form 


merical procedure. We start by determining the equilibrium 
order-parameter profile [</) * (z), 0 < z < L\, solving iteratively 
Eq. ( |4.1[ ). The solution of this equation depends, however, 
for a given range of parameters T and A/i, on the choice of 
the initial state of the order-parameter profile. The two basic 
initial states are 


z=0 


/3®ex=£) - In [l — 0 *(z) 2 ] --ln[l-( 


1 + r(z) 


+ -<b*(z)ln ~ ' T )+ --0*In 
y ' |_i — 0*(z) J v b 

+ */3AV(z)0*(z)-^[0*(z)-< 


1 + 0 * 




Pco re g, (4.4) 


i: a "liquid-like" state in which all the sites of the lattice are 
occupied by a fluid particle, i.e. [0* (z) = 1, 0 < z < ZJ; 


ii: a "gas-like" state where 0*(O) = 1, 0*(L) = 1, (j>*(z) = 
0 , 1 <z<L-l. 


where 


P OJreg — 


(a-l)K ^ 1,c+S2 ^ ~~ + a ^ Kk 


xL 


-(7+1 e o—d 


(4.5) 


As it is clear from Eq. the model presented above does 
not account for the direct wall-wall interaction, but only for 
the walls-fluid and fluid-fluid ones. Thus, in order to obtain 
the complete net force acting between the surfaces bounding 
the fluid one has to add to the force calculated from Eq. ( |4.4| > 
via Eq. the one due to the direct wall-wall interaction. 
Then the resulting net force will be the total force f lol (L\'l\ fi) 
betw e en th e pla tes b ounding the fluid. With the use of Eqs. 
(2.1 1 , (A 6 i and (4.4 1 one can write f tot (L\T,fi) in the form 


ftot(L\T,H) = [C0ex(T + l|T,jU)-C0 e x(T-l|T,jU)] 


4 Ks^ c S 2 ,c t —op O—d 
L ‘ 9 ret ) 


G(d. o)K~X 


(4.6) 


Thus one needs to calculate the profile starting from both ini¬ 
tial states. If the two final states coincide they provide the 
unique minimum of the functional - see Eq. ( |4.4| >. If they dif¬ 
fer, one has to check which one of them provides the absolute 
minimum of the grand canonical potential. The simplest way 
to clarify that question is to calculate P<X) ex via Eq. ( |4.4| >. On a 
given line in the (7) A/i) plane the values of the grand canoni¬ 
cal potential of the two profiles coincide. The set of points on 
this line defines the phase diagram of the system. An alterna¬ 
tive and practically more effective approach for determining 


of the phase diagram will be presented in Sec. VI 


Heretofore for systems governed by dispersion interactions 
there are no closed-form analytical expressions about the scal¬ 
ing function of the Casimir force even within mean-filed the¬ 
ory. Such expressions exist, however, for systems governed by 
short-range interactions GH. For (+,+) boundary conditions 
at A/i = 0 the corresponding results are 


(i) Xg as (x, > 0,x^ = 0) = - [2 K{k)f k 2 (l-k 2 ), (4.10a) 


where the last term represents the direct wall-wall interaction. 
On the other hand, if one subtracts from the potential fi) ex its 
regular part 0 J reg , i.e., if we consider the quantity 

A CO lim [(£2 £2 re g)/,0^] TcObuik, (4.7) 


with x t = [2K(k)] 2 (2k 2 — l); 

(ii) A+ r as (0 > x, > - % 2 , = 0) = -4 K 4 ( k), (4.10b) 

with x t = [2K(k)] 2 ( 2k 2 — l); 


then, in accord with Eqs. ( 2,1 3[ > - ( |2. 1 5[ ), the L dependence of 
A to via Eq. ( |2. 1 1 provides the singular part of the total force, 
i.e., the critical Casimir force - ,/cas(T|7’./i). Explicitly, one 
has 


(3E AjC as(T|r,/t) = [Aa(L+l\T,n) 
-Aa)(L-l\T,n)]. 


(4.8) 


Obviously, near T c the total and the Casimir force are related 
via the expression 


/ tot (L|r,/i) =/ 37 +Cas- 


47fsi iC S2, e 
G(d, o)K~X 


oto—d 




+ 


^-(si,c + S2,c)~]:G{d,a)Kk 

K r 4 




<7 z o—d 


(4.9) 


Eqs. (2.8 1 , (4.1 1 - ( |4.9[ > provide the basis for our numerical 
treatment of the finite-size behavior of the total as well as of 
the Casimir force in fluid systems governed by van der Waals 
interactions. Let us briefly outline the main steps of this nu- 


(iii) X£ r as (*, < = 0) = -AK\k) (1 - k 2 f , 

(4.10c) 

with x t = — [2K(k)] 2 (k 2 + 1), where K(k) is the complete 
elliptic integral of the first kind, 0 < k < 1. In Eq. ( |4.10| ) 

the scaling variable of A/ as (x t . x^ = 0) is x t = t (L/^+) = 

t (L/^q) 2 . We note that A/ as (x t , x^ = 0) is analytic for all val¬ 
ues of x t , because the film critical point [T = 7) / . p = jl, j] is 

located off coexistence at /t fj z, — jJ. c ~ l- 3 G2HT3SHIB1IIS21. 

Obviously one has x t > 0 if k > 1 / +2 (with k = 1 / +2 corre¬ 
sponding to the bulk critical point), 0 > a> > — ?r 2 if 1 /y/2 > 
k > 0 (with k = 0 corresponding to the actual film critical 
point), and x t < — n 2 if — 1 < k < 0 (negative k describes the 
region below the bulk critical point). For the field dependent 
scaling function of the Casimir force Xc as (x t ,x^ / 0) there are 
no analytical expressions even for short-range systems. 

In the next section, based on the results reported in sections 
[HI]- |IV| we present numerical results for the behaviour of the 
scaling functions of the critical Casimir force and of the van 
der Waals term for the cases d = O = 4 and d = O = 3. 
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x,=t(L/f 0 + ) 2 






FIG. 3. Behavior of the scaling function Ac r i t (x f ,xji) in a d = 3 dimensional confined fluid system. On (a) and (c) the temperature dependence 
°f Xcrit(*»!■*«) is shown at Afl = 0.0 (i.e., at x^ = 0.0), while on (b) and (d) - the field one at T = T c ,( i.e., at x, = 0.0). In (a) and (b) the 
parameters characterizing the interactions in the systems have values ,5] c = 1.0, ,? 2 ,c = 0.0 and A = 1.0, while in (c) and (d) one has ,V| c = 1.0, 
S2, c = — 1.0 and A = 1.0. The considered separations between the confining walls are L = 20 (~*~), 50 (“•"), 100 (“♦“), 200 (-*-), 400 (-T 

— ) and 800 (“★") layers. On every sub-figure the scaling functions are compared to one for a system with L = 800 layers (-) governed 

by pure short-range interactions (A = = 0.0, i = 1,2). One observes that when S2, c = 0-0 all scaling functions of systems with L > 20 are 

indistinguishable from the one for a short-range system [see (a) and (b)] and correspond to attractive force. When S 2 ,c = —1.0 the scaling 
functions do trace separate curves for different film thicknesses L. For Afl = 0.0 [see (c)] one notices that for L < 100 the scaling functions 
correspond to repulsive forces. For L = 200 the force changes sign twice. When L is further increased the force becomes entirely attractive 
for L > 800. At / — 7' [see (d)] some of the scaling functions change sign twice in the "gas" phase of the fluid medium (those for L < 100). 
Upon increasing the separation the change occurs only once, and is not observed for L > 800. When fl = fl c the curve for the short-ranged 
interactions depicts the numerical evaluation of the exact analytical result given by Eq. < [4~T0l >, while for T = T C this curve is evaluated within 
the presented mean-field theory. 


V. RESULTS FOR THE BEHAVIOR OF THE FORCES IN 
d = 3 CONFINED FLUID SYSTEM WITH <7 = 3 
DISPERSION POTENTIALS 


In the current section, using the results for d = <7 = 4 from 
the mean-field type numerical study of the system introduced 
in sections III - IV we will present some approximate results 


for the behavior A' cn |, the Hamaker term and the resulting crit¬ 
ical Casimir force, between the two surfaces confining a van 
der Waals type fluid in d = 3. 

Within the mean-field theory the dependance of the corre¬ 
sponding forces from the temperature T and the chemical po¬ 
tential difference Afl near the bulk critical point (T = T c . Afl = 
0 ) is given through the dimensionless temperature and field 


scaling variables x t = t (L/^YY and x^ = /3 c Ap (L/<jj o,q) 3 - 
In our numerical treatment we take these variables to range 
in the intervals: x, £ [— 24 2 ;24 2 ] = [—576;576] and x^ € 


[—24 3 ;24 3 ] = [—13824; 13824]. We consider systems with 
fluid layers L being 20, 50, 100, 200, 400 and 800. We fix 
one of the walls-fluid coupling parameters (say ,V| j to have 
a value either 1.0 or —1.0, while the other one, S 2 , c , is varied 
from 0.0 to — 1.0 with a step of —0.2. The fluid-fluid coupling 
parameter A is supposed to be either 1.0 or 2.0. 

In order to determine the behaviour of the above mentioned 
quantities in a d = 3 confined fluid system, using the presented 
approximation, we proceed in the following way. First, within 
the mean-field treatment of the critical behaviour {d = <7 = 4) 
we solve the equation for the equilibrium order parameter pro¬ 
file - Eq. ( |4.1| >. Next, with the use of Eqs. ( |4.4| i and ( |4.6[ > we 
obtain the total force of interaction between the plates. Sub¬ 
tracting from Eq. ( |4.6| > the Hamaker term Eq. (|A7[ ) we end up 
with the scaling function Y cr j t [see also Eq. ( ]2.8 [ >] . In order 
that this function contributes properly to the critical Casimir 
and hence to the total force of interaction in d < 4, one must 
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FIG. 4. Behavior of the scaling function X c rit(xr,Jqi) in a d = 3 dimensional confined fluid system. On (a) and (c) the temperature dependence 
°f Xait( x t,Xu) is shown at Ap = 0.0, (i.e., x^ = 0.0), while on (b) and (d) - the field one at T = T c , (i.e., at x, = 0.0). In (a) and (b) 
the parameters characterizing the interactions in the systems have values i'i c = —1.0, S 2 ,c = 0.0 and X = 1.0, while in (c) and (d) one has 
.s | c = —1.0, S 2 ,c = —1.0 and X = 1.0. The behavior of the scaling functions in a system characterized by X = 1.0, ,S| c = —1.0 and S2,c = 0-0 
is similar to that of such with X = 1.0, si jC = 1.0 and ys.c = —1.0 [compare wit Fig. |3] (c) and (d)]. For X = 1.0 and si jC = S 2 . c = —1-0 at 
Ap = 0.0 [see (c)] one observes that for L = 20 (-•-) the scaling function changes sign twice, having two minima and a maximum. When 
the separation L is increased the values of the minima decrease rapidly towards zero, while that of the maximum increases, being highest for 
L = 100 (“♦“). When L = 200 (-*-) the scaling function corresponds to repulsive force, now having two maxima and a single minimum. 
Reaching L = 400 (~T~) the scaling function changes sign twice, and corresponds to attractive force forL > 800 (“★-)• At T = T C [see (d)] for 
L = 20 and 50 the scaling functions change sign once for x^ > 0 [see (d / )]. Then for L = 100 one observes triple sign change, with a shallow 
minimum in the region x^ > 0 and pronounced one in x^ < 0. For L = 100 the maximum of the scaling function has its highest value. Upon 
increasing the separation sign change of the scaling function occurs twice, and is not observed for L > 800. Note that for L = 800 the scaling 
function still changes sign in the region x^ < 0. When p = fl c the curve for the short-ranged interactions depicts the numerical evaluation of 
the exact analytical result given by Eq. (|4.10[>, while for T = T C this curve is evaluated within the presented mean-field theory. 


normalize it accordingly. The need of such normalization is 
explained in details in Ref. l76l (see there Sections IV.A. 1 
and IV. A.3). For boundary conditions T one has 


* 2 <-)= 


2AW(d = 3) 

*2;r('=AM = o) 


« 0 + ( 0 ) 


n 4 


(A)J 


*irit’ MF (-), (5-1) 


where is the value of the scaling function for a system 

within mean-field treatment governed by short-range interac¬ 
tions at its corresponding bulk critical point, and X ;MF (.) 
is the scaling function of the critical Casimir force, calculated 
for d = <J = 4, with A =/=■ 0, si iC 0, S 2 ,c / 0- Here A( T ) (d = 3) 
is the Casimir amplitude for the d = 3 Ising universality 
class with boundary conditions t, while (0) and (A) 
are the amplitudes of the bulk correlation length in mean- 


field systems with, correspondingly, short-ranged (A = 0) and 
long-ranged (A 0) fluid-fluid interactions, as (A) = sjvi 
(see Eqs. (4.15) and (4.17) in Ref. l76l ). Therefore for 
A = 0, 1 and 2 one has t ; 2 = 1/9, 0.1640, and 0.1998, and 
hence £q" = 1/3, 0.4050, and 0.4470 respectively. Taking 
into account that the value of the Casimir amplitude for the 
d = 3 Ising universality class with (+,+) boundary condi¬ 
tions □sa is 


A (+ ’ +) = -0.410(29) (5.2) 

and that X^/^ F (r = Aji = 0) = —1.7315, for the normalizing 
coefficients one obtains: 0.217 for A = 1 and 0.147 when A = 
2 . 


The dependence of the scaling function A cnl on the temper¬ 
ature x t and field scaling variables is summarized in Figs. 
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FIG. 5. Interplay between the scaling function X c rii( x t> x n) (-) and the van der Waals term in excess of its regular contributions 

2/3 H& ex (x t ,Xn) (“♦“), resulting in the scaling function of the critical Casimir force /3F)i,Cas(*f Afi)^ 3 = Zc as — ^cnt(xt,x^) + 2pH^ ex (x t ,Xn) 
(-+-) in a d = 3 dimensional confined fluid system at A fl = 0.0, (x^ = 0.0) [(a) and (c)] and T = T c , (xt = 0.0) [(b) and (d)]. The separa¬ 
tion between the walls confining the fluid medium is taken L = 20 and the parameters characterizing the interactions in the systems have the 
following values: A = 2.0, S2, c = 0-0; [(a) and (b)] = 1.0; [(c) and (d)] = —1.0. 


[3] and [4] for different separations L. 

Once having a good approximation of Z cr i t in d = 3, ob¬ 
tained in the way described above, one can determine the scal¬ 
ing function of the critical Casimir force of the fluid system by 
adding to Z cr it the singular part of the Hamaker term given by 
Eq. ( |A8| >, i.e., that part of it that depends on </>/,. When per¬ 
forming this procedure one must pay attention to the use of 
the coupling parameters s, jC , i = 1,2, because they also de¬ 
pend on the dimensionality d of the system. Using Eqs. (3.31, 
( |3.7| > and ( |3. 11 [ i one has 


,, G(3,3) 

s Uc {d- a- 3)— G(4 , 4 )S«, 


(d = cr = 4), i= 1,2. (5.3) 


As it is clear from the above, in system governed by short- 
range interactions X cnt coincides with Z f s .' as is negative, 
which corresponds to attractive critical Casimir force, for any 
value of the scaling variables x, and under (+,+) bound¬ 
ary conditions - see Eq. ( |4. 1 ()[ > for X^ as in the case x^ = 0, as 
well as Fig. 0 a) and (b). Note that under such boundary con¬ 
ditions one observes an enhancement of the order parameter 
due to the confinement in comparison to the order parameter 
at the same distance from an individual wall. On that ground 
in systems with long-ranged interactions such that ,V| ,c>0 and 
—,V| . ( <C S 2 .c < 0 we expect X CI n to remain again negative for 


any separation L and at any value of the scaling variables, ir¬ 
respective of the value of A (which is always non-negative). 
Indeed, this turns out to be true and is depicted in Fig. 0 a) 
and (b). However when S 2 . c ~ — .V| x and the separation be¬ 
tween the walls is relatively small a significant part of the 
system is disordered which results in non-negative or sign¬ 
changing scaling function [see Fig. 0c) and (d)]. As the 
distance L is increased the influence of the effective surface 
potential AV(z) [see Eq. ( |3.6[ ) | quickly decreases and only the 
additional ordering effect of the fluid-fluid interactions influ¬ 
ences the behavior of the order parameter and hence of Z cr j t . 
Along with the same line of arguing, when both wall-fluid 
coupling parameters are negative, Z ml is negative for any x t 
and Xfx only for very large separations L where the effect of 
the long-ranged interactions on the behavior of the system is 
negligible. Naturally, since the short-ranged surface poten¬ 
tials do support (+,+) boundary conditions, the role of the 
negative substrate potentials, which oppose the order near the 
boundary, will be stronger than that of the positive substrate 
potentials which try to reinforce the phase preferred near the 
boundary. For example, we observe that the behavior of Z cm 
in a system with si c = 1.0, S2,c = —1-0, A = 1.0 and in such 
with si c = 0.0, S 2 , c = —1.0, A = 1.0 is almost identical for 
any L, both as a function of x, and x^ [compare Fig. 0c) 
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and (d) and Fig. [dj a) and (b)] 11611 . Thus, for a fixed A 
the behavior of the scaling function is mainly determined by 
the interplay between the short-ranged surface fields and the 
strong negative wall-fluid coupling « 2 ,c- If H,c = $2.c = — 1.0 
and A < 2.0 the scaling function A ml exhibits an unexpected 
behavior as a function of L: for moderate values of L the max¬ 
imum of the repulsive part of the force increases with increas¬ 
ing L both as a function of x t and [see Fig. |4jc) and (d)]; for 
larger values of L the maximum decreases, as expected, and 
the overall behavior of the scaling function approaches that 
one of the system with completely short-ranged interactions. 

The behavior of the critical Casimir force Aq is is depicted 
on Fig. [5] In order to illustrate only the main idea, here we 
restrict ourself to the choice of parameters A = 2 . 0 , S2, c = 0-0 
while ,v 1 c = ±1.0. The data for the Casimir force and the 
Hamaker term are presented for a system with L = 20 layers. 
First, let us recall that for T >T C at A fi = 0 one has 0* = 0 and, 
hence, the behavior of Xcas and that of coincide in this re¬ 
gion [see Fig. [5}a) and (c)]. On the other hand, if T <T C one 
has (j>b 7 ^ 0 , with the singular (i.e., the (pi, dependent) part of 
the Hamaker term corresponding to repulsion for ,V| ,C> 0 and 
attraction otherwise, see Eq. ( |A 8 | ). For ,v 1 c = 1.0 the scaling 
function A cr i t is negative, but the singular part of the Hamaker 
term is positive, i.e., repulsive in the low-temperature region. 
Thus, the Casimir force is attractive near and above T c but be¬ 
comes repulsive below T c . For si c = —1.0 one has that is 
non-negative, but with the singular part of the Hamaker term 
being negative, i.e., attractive for T<T C . Thus, the resulting 
critical Casimir force Yeas changes sign from being repulsive 
above T c to becoming attractive slightly below T c . The be¬ 
haviour of the Casimir force as a function of x^ at T = T c 
for the two sub-cases of si c = ±1.0 is depicted on Fig. |5jb) 
and (d). For ,v 1 c =1.0 one observes that A cnt < 0, while the 
singular part of the Hamaker term changes sign from positive 
to negative with All decreasing. Thus, the resulting Casimir 
force changes sign once from being slightly positive (repul¬ 
sive) for All > 0 to negative (attractive) for All < 0. The case 
ii jC = — 1.0 is much more interesting. Since then A cr j t changes 
sign twice for All < 0, while the singular part of the Hamaker 
term changes sign ones, the resulting Casimir force happens 
to change its sign tree times as a function of x u . 

One might pose the question about the phase behavior of the 
thin fluid systems with such competing surface and substrate 
potentials. We consider that question in the next section. 


VI. PHASE BEHAVIOUR OF SYSTEMS WITH 
DISPERSION FORCES 

The phase behaviour of a confined fluid medium between 
parallel walls exerting identical surface adsorption poten¬ 
tials i)/i| = y s V J = Sfj .2 = Jsr’ 1 7 ^ 0 on both surface lay¬ 
ers has been studied extensively both theoretically and ex¬ 
perimentally ESI cm ESI Ell E3 Ell ESI E3 H371 
1162111631 . So far the considerations have been made either 
for fluid systems governed by pure short-range interactions 
01341114311144111461114711157111631 . or for systems with sur¬ 
face potentials strongly preferring the same phase of the fluid 
ll76l [1361II371II621 . Here we are going to consider the case 
when at least one of the wall-fluid potentials favor a phase of 
the fluid different from the one preferred by the short-ranged 



FIG. 6. Phase diagrams of a d = 4 confined finite-size fluid system 
with strong adsorption of the "liquid" phase of the fluid medium on 
the confining walls, for different separations L between them. The 
values of the coupling parameters characterizing the interactions in 
the system are = 1.0, 57, c = —1-0 and A = 2.0, and the coordi¬ 
nates of the finite-size critical points for the different separations are: 
bottom-up L = 20 (▼) (0.98779783,-1.8421 x 1(T 3 ); L = 30 (A) 
(0.99553852,-0.5957 x 10~ 3 );L = 40(*) (0.99743470, -0.2684 x 
1(T 3 ). 


surface potentials. 

The phase behavior of systems with L = 20,30,40 and of 
the bulk system with si jC = 1.0, 52 ,c = —10 and A = 2.0, ob¬ 
tained within our mean-field model, is presented in Fig. [6] 
A relatively detailed study of the influence of the different 
values of si c , S 2 , c and A on the phase behavior of thin fluid 
films (in the case considered L = 20) is illustrated in Figs. 
[7] and [9j where Fig. [9] represents blow-up views of some of 
the phase diagrams very close to the capillary condensation 
critical points T c l . We observe that the overall phase behav¬ 
ior of the system is similar to that one obtained for a short- 
ranged system by Fisher and Nakanishi II 1441 - one has a line 
of a capillary condensation shifted towards the gas bulk phase 
which terminates at its own finite-size critical point T c p. One 
can identify the following properties of the phase diagrams 
presented in Figs. [ 6 ][7] an 

i) The precise position of the phase diagram depends on L. 

This property is well known also for systems with short- 

ranged interaction. 

ii) For a fixed L the precise position of the phase diagram 
depends on ,v | c , S 2 , c and A - smaller L stronger the corre¬ 
sponding influence of these parameters. 

Negative wall-fluid couplings .v 1 c and/or S 2 . c have disorder¬ 
ing effect on the system which manifests in increase of the gas 
phase and decrease of the liquid one near the surfaces. Gen¬ 
erally the phase diagram, as then expected, appears above the 
one of a short-range systems [see Fig. [7jb)] since in order to 
"liquify" the system one now needs higher "pressure". 

iii) The position of 7) / , thus the shift of the critical point 
of the finite system with respect to that one of the bulk one, 
depends not only on L, as in the short-ranged systems, but also 
on .?i c , «2,c an d A. 

For A 7 ^ 0 the ordering in the system increases and, hence, 
one observes increase in both the critical temperature shifts 
as well as in the chemical potential difference compared to 
systems governed by pure short-range interactions. When 
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FIG. 7. Comparison between the phase diagrams of d = 4 confined 
finite-size fluid systems with strong adsorption of the "liquid" phase 
of the fluid medium on the confining walls, obtained for L = 20 
and various values of the coupling parameters: (a) bottom-up (A) 
s\ c = 1.0, S 2 ,c = 0.0, X = 2.0; (■) c = 1.0, 52 ,c = 0-0, X = 1.0; 
(•) Si.c = s 2 ,c = x = 0.0; (T) si, c = 1.0, 52, c = -l.o, X = 2.0; 
(^) *t,c = 1-0, 52, c = —1-0, X = 1.0; (b) bottom-up (•) 
Ji,c = s 2,c = X = 0.0; (X) 5i, f = -1.0, 5 2 , c = 0.0, X = 2.0; (V) 
5i,c = -10, 5 2 , c = 0.0, X = 1.0; (<))5 l c = -1.0, 5 2 , c = -1.0, X = 
2.0; (4k) 5i, c = —1.0, 52 , c = —1.0, X = 1.0. The coordinates of the 
finite-size critical points [T c ^/T c , — fJ-c) /(7b7c,l)] of the differ¬ 
ent systems are: (a) bottom-up (A) (0.98912146,-2.6691 x 
10 3 ); (■) (0.99046179,-2.1410 x 1(T 3 ); (T) 

(0.98779783,-1.8421 x 10- 3 ); (•) (0.99376528,-1.1101 x 10- 3 ); 
(♦) (0.99185806,-1.2482 x 10 3 ); (b) bottom-up (4) 

(0.99130666,-1.5386X 10~ 3 ); (•) (0.99376528,-1.1101 x 10 3 ); 
(^) (0.99279091,-1.0222 x 10~ 3 ); (0) (0.99303773,-0.9915 x 
10~ 3 ); (4b) (0.99468863,-0.5085 x 10~ 3 ). The phase diagram of a 
bulk system (L —> 00 ) is shown as thick black line ending at the bulk 
critical point (•). 


and 52 , c are both positive this effect is even further enhanced. 
If, on the other hand, for fixed X ^ 0 the parameters 5i, c and 
52 , c are both negative or, at least, have different signs, the 
shift of the finite-size critical point with respect to that one in 
short-range systems can be completely canceled or, for strong 
enough potentials, this shift can even further decrease via T c l 
moving toward the corresponding bulk critical point. 

iv) The competing effect of T, A fi, 5i c , 52 , c and A on the 
ordering in the system leads near the capillary condensation 
critical point for thin films to a non-monotonic behavior of 
the phase diagram - see Fig. [9] The behavior of the system 



FIG. 8. The scaling functions of the total isothermal compressibility 
KjLT^I v [(^), (A) and (▼)] and the phase diagram (■) of a d = 4 
finite-size fluid system, confined between parallel walls both strongly 
adsorbing the "liquid" phase of the fluid medium, for L = 20 and 
values of the coupling parameters 51 c = 1.0, 52 , c = 0.0 and X = 
1.0. The bulk critical point is denoted by the symbol •, while the 
finite-size critical point (T = T c ^, fl = jl c ,i) is marked by the symbol 
■ . The curve designated by (^) shows the behavior of Kj-L~7/ v at 
jit = /If. This behavior at the finite-size critical point (T = T c i,fl = 
LI c ,l) is depicted with the symbol (A). When T < T CiL and fl < 
Ll C L, the compressibility Kj changes with a finite jump (▼), being 
an indicator of a first-order phase transition. The points of first-order 
phase transition form the coexistence curve, i.e., determine the phase 
diagram in the phase plane. 


with at least one of the wall-fluid potentials being negative is 
especially interesting - see Fig. |9]d), where in a small tem¬ 
perature interval at phase coexistence when T decreases Att 
increases. 

Below we briefly explain how the above phase diagrams 
have been calculated. 

On Fig. [ 8 ]we illustrate the method used to obtain the phase 
diagram and the finite-size critical point at which the coex¬ 
istence curve ends of finite systems characterized by long- 
ranged fluid-wall potentials. We do that by studying within 
our mean-field model the scaling functions of the total isother¬ 
mal compressibility KjLT^I v . When /i = jtt c one has that 
k'/ L F 1 ' is a smooth function of the temperature, even in the 
vicinity of the bulk critical temperature, reaching its maxi¬ 
mum for T >T C . For < /i < ji, the value of the maximum 
increases, being still finite, and appearing at T < T c . Upon 
reaching the finite-size critical point ( T = 7, / . ,u = /I r /J, 
K T L-r/ v diverges, giving rise to power-law singularity. At 
temperatures T < T r j and chemical potentials jx < /i r / , the 
value of KjLT^! v changes with a finite jump, being an in¬ 
dicator of a first-order phase transition. The points of first- 
order phase transition form a line - the coexistence curve in 
the phase plane, separating the observed phases: spin-"up" 
phase in magnetic systems, liquid phase in one-component 
fluid systems, A-rich phase in binary liquid mixtures and, cor¬ 
respondingly, spin-"down" phase, gas phase, B-rich phase. At 
the critical point and beyond the differences between these 
phases disappear. The resulting phase diagrams are shown in 
Fig. □ where (a) summarizes the curves with 51 = 1.0, and 
(b) - those with 51 = — 1 . 0 , where different combinations of 
values of 52 and A are considered. The line marked with filled 
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FIG. 9. Blow-up views of some of the phase diagrams at d = 4 very close to the capillary condensation critical point T c i for L = 20. The 
influences of the different dispersion interactions (walls-fluid and fluid-fluid) is visualized. The phase diagram of a short-range system is 
given in (a) in order to serve as a reference on which the changes due to the influence of fluid-fluid interaction, i.e., on A - see (b), or on 
the wall-fluid interaction, i.e., on s, - see (c) and (d) is studied. In (b) only the influence of the long-range fluid-fluid interactions is shown 
(A 7 ^ 0 , .si jC = ^ 2 ,c = 0 . 0 ), while on (c) and (d), that of a single wall-fluid is investigated ( 5 i jC ^ 0 , A = 52 , c = 0.0). 


(red) circles represents the phase diagram for a system with 
completely short-ranged interaction. 

As it was shown in Refs. 1 1361 1371, for the total isothermal 
compressibility per particle (or susceptibility, if one considers 
magnetic systems), one has 


k t 


1 

L + l 


I>r(z) 


L+l 



( 6 . 1 ) 


where R 1 is the inverse of the matrix R with elements 


R 


1 - 0 2 (z') 




(6.2) 


and 


respect to the field h(z*) = /3 [A/r — AV (z*)]/2. 

In order to determine Kt(x,,Xij \L, {s,-, i = 1,2}, A) or its 
"scaling function" X K = L ' ! + xy in a fluid film with thickness 
L, we first solve numerically Eq. ( |4. 1 [ 1 within our mean-field 
model, which allows us to determine the matrix R with the 
use of Eq. d^2l >. After that, using Eq. s we obtain the 
"local" isothermal compressibility Ky (z), and summing over 
z, the total isothermal compressibility. We recall that within 
the mean-field treatment one has v = 1/2 and 7=1. 

The phase diagrams depicted on Fig. [TJresult from the inter¬ 
play between the fluid-fluid and walls-fluid interactions. Note 
that the increase of A at fixed s LC . i = 1,2 lowers p A/i at co¬ 
existence, while the decrease of the wall-fluid coupling at one 
of the walls (say S 2 . c ) keeping A and ,v 1 c constants, increases 
it. 


K- r (z) = £G( Zj z*) = £( R - 


' Z,Z 


= 4 I^(0,^(r|,z*))-(^(0,z))(^(r|,z*)) 


(6.3) 


VII. EXPERIMENTAL FEASIBILITY OF THE 
PREDICTED EFFECT 


is the "local" isothermal compressibility, which reflects the 
response of the system from a given layer due to the change 
of the external field in that layer. More precisely, in Eq. 
G(z,z*) is the density-density correlation function, G(z,z*) = 
8(j)(z)/[28h{z*)], where the functional derivative is taken with 


In the fast emerging field of micro- and nano- devices, the 
fluctuation induced and dispersion forces like the Casimir and 
van der Waals ones, associated with the interaction between 
the individual components of such devices, play an essential 
role. Due to the quantum origin of these forces, currently there 
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are no theoretical as well as experimentally suggested control 
parameters that can be used to conveniently and reversibly 
govern the sign and the strength of these forces if the inter¬ 
acting objects are in vacuum. The situation changes if these 
interactions take place in a fluid. Very recently in Ref. HSU 
the authors suggested controlled quantum Casimir levitation 
caused by the introduction of thin film surface coating on a 
porous planar substrates both immersed in a suitable fluid. 
More specifically, they have shown that the Casimir force be¬ 
tween Teflon and cassiterite (SnC>2) nanosheet immersed in 
cyclodecane is attractive at large separations but repulsive at 
small separations, resulting in a stable equilibrium distance 
where the total force is zero. In Ref. ED the authors re¬ 
ported the use of aerogels, yielding repulsion down to sub¬ 
micron distances at realistic porosities. The described so far 
set-ups are quite close in spirit to the considerations presented 
in the current article with the basic difference that the force 
will be due to the critical fluctuations of the fluid, instead of 
the fluctuation of the electromagnetic field. We suppose that 
the usage of low density substrates like aerogels can be also 
applied in the thermodynamic Casimir effect since the param¬ 
eters Sj directly depend, as shown in Eqs. ( |3.7[ ) and ( |3. 11 [ i, on 
the densities of the slabs. 


TABLE I. Physical characteristics of the considered confined fluids 
(column 1) - the critical density p c in units g/cm 3 (column 2), tem¬ 
perature T c measured in K (column 3), pressure P c in MPa (column 
4), the lattice constant oq in A (column 5), the distance ri also in A 
(column 6) at which the inter-particle potential is zero and the poten¬ 
tial well depth J l st in units k^T c (column 7). The values of p c , T c and 
P c for argon are taken from Ref. (165) . for krypton from Ref. 11661 
and for xenon from Ref. (163. The value of qq is calculated based 
on the data for p c (see the text) and the presented values of rL and jj. r 
are taken from Ref. [ 1681 . 


1 
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4 

5 

6 

7 

fluid 

Pc 

T c 

Pc 

ao 

J 
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Ar 

0.536 

150.663 

4.86 

4.98 

3.867 

0.798 

Kr 

0.908 

209.480 

5.53 

5.35 

4.165 

0.790 

Xe 

1.113 

289.765 

5.84 

5.81 

4.512 

0.794 


In order to show with concrete examples that the exper¬ 
imental observation of the effects theoretically predicted in 
the current article is possible, we consider the occurrence of 
critical fluctuations in simple non-polar fluids such as argon 
(Ar), krypton (Kr) or xenon (Xe), confined between graphene 
coated carbon nanotube aerogel substrate (CAS) ED and a 
silver (Ag) or gold (Au) one. Both confining surfaces are as¬ 
sumed coated by monatomic (or very thin) lead (Pb) film, en¬ 
suring the (+,+) boundary conditions. The basic data needed 
for the calculations concerning the considered fluids are pre¬ 
sented in Table [T] Furthermore, from the data reported in 
Ref. 068] we have that /3 c J s f’ Pb = 1.79, /3 c /£ r,Pb = 1.32 and 
/3 c J* e,Pb = 0.96, where = 1 / (IcbT c ), see Table]]] is specific 
for each fluid medium. As far as the geometry of the experi¬ 
mental set-up matters, one can choose to study the predicted 
behavior of the force either in plane parallel geometry Es¬ 
ina, or consider such where one of the plates is flat and the 
other has the shape of a spherical lens with large (> 1 cm) 
curvature radius 11731 . As pointed out in Ref. 11711 . the diffi¬ 


culties in the usage of the plane parallel geometry are bigger, 
which leads to lower accuracy of the results in comparison 
to the other configuration. The use of plates in the form of 
spherical lens to study effects predicted for systems composed 
out of flat parallel surfaces can be mathematically justified, as 
usual, by the Derjaguin (151 ;, T741 or the gradient expansion 
1175441771 approaches, which are applicable for calculations 
of the interactions between curved objects. 

In what follows we evaluate the coupling parameters s c 
and A which are the basic parameters reflecting the mate¬ 
rial properties in the theory for the behavior of the force pre¬ 
sented above - see Figs. [3] — [5] From Eqs. ( |3.7[ ) and ( |3.11| ) 
we have that s c = 0.5G(d,o)[p n dp c J'' s — p c j3 c /'], and thus, 
in order to determine s c , one needs to know the long-range 
inter-particle interaction energies J Lx and J l , as well as the 
number density p m j (see column 7 in Table |II} of the confin¬ 
ing substrate relative to that of the fluid medium p c . Since 
within the mean-field theory the number density of the fluid 
is p c = 1 /2 one has p m j = 0.5p[mol/cm 3 ]/p c [mol/cm 3 ], with 
p [mol/cm 3 ] = p[g i /cm 3 ]/ (A T uN a), where A r is the standard 
atomic weight of the substance, u - the atomic mass unit and 
Na is the Avogadro constant. From here one can also estimate 
the lattice constant of the fluid medium at the critical point; 
one has ao [A] = 10 10 x {p [g/cm 3 ]/(A r M)} -1 / 3 (see column 5 
in Table |TI]). We will make the assumption that all interactions 
between the constituents of the system are of Lennard-Jones 
type, i.e., the interaction potential can be written in the form 



(7.1) 


where r m = 2 1 / 6 ro is the distance at which the potential 
reaches its minimum, ro is such that vvu(ro) = 0 and ,/ sr is the 
depth of the potential well. Therefore one has that J LS = 2J sr 
(see columns 5 in Table |fi|) and A = J 1 /J l ST = 2 (see column 7 
in Table|ll. We evaluate the cross interaction potential param- 

l l S 

eters 7 sr and r 0 ’ based on Kong’s combining rules, because of 
their appreciated accuracy II1781 



(7.2a) 


jl,S 
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Wo) 

2 13 


12 


1 + 
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L4K) I2 J 


121 1 / 13 ' 


13 


(7.2b) 


After all that we are ready to compute the key parameter s c 
that influences the behavior of the force for moderate values 
of L. The obtained results are reported in column 8 of Table 
E • One observes that, depending on the choice of materials, 
one can indeed have both positive and negative parameters 
s c . The existence of combination of materials for which s c is 
negative represents the main prerequisite for the experimental 
relevance of the effects predicted in the current article. 

Although the mean-field theory gives poor quantitative es¬ 
timation of the behavior of the scaling function Xc as , it is 
tempting to, nevertheless, evaluate the corresponding force 
/pas based on it. To do so we consider the interaction be- 
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TABLE II. Physical characteristics of the confining silver (Ag), gold (Au) and carbon based (C) substances and of their interactions with the 
fluid media - the distances and in A (columns 2 and 4) at which the inter-particle potential within the substrate and between it and the 

fluid is zero, the corresponding potential well depths 7j r and j' s f in units k^Tc (columns 3 and 5), the density p of the substrates in g/cnr 3 
(column 6), the number density p n( j (column 7), the substrate-fluid coupling parameter evaluated at the critical temperature - s c (column 8 ) 
and the Young modulus E of the corresponding substrate in GPa (column 9). The values of r q and 7* r are taken from Ref. fl68l . while those 
of rij s and J l s f are calculated via Eqs. <7.21. The density of the CAS (the bottom value in column 6) as well as its Young modulus (the bottom 
value in column 9) are taken from Ref, fl 691 . 


1 

2 

3 

4 

5 

6 

7 

8 

9 

fluid/substrate 

r s 

r 0 

PcJl r 

Ls 

r o 

Pc4? 

P 

Pnd 

«c 

E 

Ar/Ag 


0.120 

3.597 

0.258 


3.62 

1.12 


Kr/Ag 

3.148 

0.087 

3.804 

0.194 

10.49 

4.49 

1.00 

83 

Xe/Ag 


0.063 

4.061 

0.142 


5.74 

0.89 


Ar/Au 


0.130 

3.628 

0.278 


3.65 

1.29 


Kr/Au 

3.239 

0.094 

3.831 

0.211 

19.30 

4.52 

1.17 

79 

Xe/Au 


0.068 

4.083 

0.156 


5.78 

1.06 


Ar/C 


0.351 

3.864 

0.525 


0.043 

-0.788 


Kr/C 

3.851 

0.252 

4.033 

0.428 

0.014 

0.054 

-0.779 

3.86 

Xe/C 


0.182 

4.245 

0.341 


0.069 

-0.782 



tween a CAS planar substrate and a gold disk with radius 
R = 4 ll m and thickness h = 2 um, immersed in xenon. We 
will assume that the separation between the confining sur¬ 
faces is L = 20ao — 12 nm. Then for ju = fJ. c one shall ob¬ 
serve double sign change of the force with the occurrence 
of two maxima, at which the force is repulsive and a mini¬ 
mum, at which it is attractive. The magnitude of the force 
at the maxima is predicted to be Teas, max! — 16 nN and 
Ecus. max 2 — 7.6 nN, observed at T ~ 283.5 K and T ~ 300 K 
respectively. The minimum is found at T ~ 291 K and its 
value is Tc as , min — —24 nN. Finally we predict that the two 
sign changes occur at T ~ 288.9 K and T ~ 295 K respec¬ 
tively. Note that the predicted values of Tc as are significantly 
larger than the weight of the disk, which is approximately 
2 x 10 2 nN. We stress that the predicted magnitude of the 
force is well in the range of the one reported in Ref. 11791 . 
where the authors estimated that an object with an effective 
interaction area szi = 1 lim 2 at a distance L = 10 nm to a wall 
experiences a force Teas — 4 nN when immersed in a fluid with 
T c = 300 K (about the critical one of a water-2,6-lutidine mix¬ 
ture). Let us also mention that although the presented study 
is limited to the interaction between parallel planar surfaces 
through some fluid medium, it can be straightforwardly ex¬ 
tended to objects of various geometries, e.g., via using say the 
Derjaguin or the gradient expansion approaches to calculate 
the interactions between curved objects. Then one can study 
also the behavior of colloidal particles immersed in a critical 
fluid or in a binary mixture close to its consulate point. Thus, 
another possible application of our findings is for governing 
the behavior of colloids in a critical solvent - imagine coated 
colloid particles. Then, under proper choice of the core mate¬ 
rial of the colloid and of the coating the particles will have a 
stable equilibrium distance from the bounding plate that will 
depend on the temperature and the pressure of the fluid. Thus, 
we hope that our result can be used as a guide to the design of 
new experiments on colloid systems. 


VIII. SUMMARY AND CONCLUDING REMARKS 


In the current article, in view of future experiments explor¬ 
ing and potential devices utilizing the critical Casimir forces 
in classical one- or two-component fluids confined by paral¬ 
lel substrates at a distance L from each other, we studied the 
behavior of the critical and the near critical Casimir force and 
its interplay with the van der Waals force, as well as the net 
force due to both of them. In the envisaged set-up the walls of 
the slabs are considered coated by thin layers exerting strong 
preference to a given phase of the fluid modeled by strong ad¬ 
sorbing local surface potentials, while the slabs, on the other 
hand, influence the fluid by long-range dispersion potentials 
at least one of which supports the opposite phase of the fluid. 
The fluid-fluid interactions are assumed to decay as r id ' °- 1 , 
with r being the distance between the particles of the medium, 
while the walls-fluid interactions decay as z~ a where z is the 
fluid particle-single wall separation. The strengths of the dis¬ 
persion interactions in the system are depicted via the dimen¬ 
sionless coupling parameters - A and s,-, i = 1,2, account¬ 
ing for the fluid-fluid and slabs-fluid interactions, respectively. 
While A is always non-negative, for the slabs-fluid coupling 
parameters si and sj we assumed that ,v i VJ 0 while sj < 0 . 
The confined fluid medium is modeled within the framework 
of the lattice gas model defined in Sec. m After developing 
a theory for arbitrary dimension d and long-range decay ex¬ 
ponent a of the dispersion interactions - see Appendix |A| and 


IV 


we 


studied in Sec. [v]the forces in d = 3 systems, 


assum¬ 


ing that the dispersion interactions are of "genuine" van der 
Waals type, i.e., d = <J = 3. Away from the critical region 
the total force / tot between the plates of the bounded fluid is 
simply proportional to the Hamaker constant - see Eq. ( |2.2[ ) 
and Eq. ( ]2.3[ i, while near (T = T c . All = 0) one has additional 
contribution to the force f tot which is proportional to L^ d X CI j t , 
where Y cn i is a scaling function reflecting the role of the criti¬ 
cal fluctuations of the order parameter - see Eq. ( |2.6| ). 

Using general scaling arguments, as well as explicit mean- 
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field model calculations, we have obtained the following main 
results: 

(1) In terms of some critical thickness L cr ; t of the thin films 
- see Eqs. ( |2.10 i and ( |2. 11 | i, we have established a criterion 
for the importance of dispersion forces within the critical re¬ 
gion of the system. For L < L cr j t the contributions of the dis¬ 
persion forces are important and cannot be neglected neither 
within the critical region of the system, including the bulk crit¬ 
ical point, nor outside of that region - see Figs. 00and[4] 
where the validity of these statements is clearly visualized. 
In the opposite case, i.e., when the separation L between the 
walls is much larger than L cr ; t the dispersion forces provide 
only corrections to the leading behavior of the critical Casimir 
and the total force within the critical region of the system. 
Outside the critical region, however, i.e., for T ^ T c . Ap ^ 0, 
the influence of the dispersion (van der Waals) interactions 
becomes essential for the behavior of the total force between 
the bounding surfaces. This is of experimental importance be¬ 
cause it is difficult to thermodynamically position the system 
right at the critical point. 

(2) In Appendix |A] we have derived expressions for the be¬ 
havior of the Hamaker term in the total force between the 
plates - see Eq. ( |3. 16| > and Eq. ( | A6 [ > . that allow for a study 
of the temperature and the field dependencies of the Hamaker 
constant including that ones within the critical region of the 
system. The behavior of this part of the total force as a func¬ 
tion of the scaling variables x, and for different values of 
ii, 52 and A is shown in Fig. [5] 

(3) We demonstrated that for a suitable set of slabs-fluid 
and fluid-fluid coupling parameters the competition between 
the effects due to the coatings and the slabs can result in sign 
change of - see Figs. [3] and [4] as well as of the Casimir 
force - see Fig. [5] when one changes the temperature T, the 
chemical potential of the fluid u, or L. The last can be used in 
designing nano-devices and for governing behavior of objects 
at small, below micrometer, distances. 

(4) In addition to the critical behavior of the considered 
forces, in Sec. [VT]we discussed and gave detailed explana¬ 
tion for the near critical behavior of the phase diagram, see 
Figs. @0 and [9] of finite-size systems governed by dispersion 
interactions. The influence of these interactions on the shift of 
the finite-size critical point has also been commented. 

(5) In the final Sec. VII in scope of possible experimen¬ 
tal realization of the theoretical set-up, we explored the usage 
of concrete substances. For each pair of substrate and fluid 
we calculated the coupling parameters s c and A, and used 
this knowledge to make estimations for the critical Casimir 
force Teas at several temperatures where it exhibits extrema or 
changes its sign. We have demonstrated that by proper choice 
of the fluid and substrate materials one can indeed achieve 
5 e > 0 and s c < 0, which is the main prerequisite for the ex¬ 
perimental feasibility of the effects predicted in the current 
article. 

Fet us stress that the problem of quantitative description of 
the mutual influence of the fluctuation of the electromagnetic 
field and the order parameter fluctuations of a medium when 
it is close to its critical point is extremely complicated and 
currently there is no general theory available to scope with 
it. The Fifshitz theory - which is the basic one for study¬ 
ing the QED Casimir effect has never been meant to nor can 
deal with the problem of a critical medium between two other 


media. The main quantity which knowledge is required for 
practical applications of this theory is the dielectric permittiv¬ 
ity e(a>). It is normally tabulated at some temperature, usu¬ 
ally the room temperature. It shall be noted, however, that in 
a critical fluid s(co) is itself a singular function of the tem¬ 
perature II181 182 1. We are not aware even of a theory that 
can reliably predict quantitatively how e(co) will depend on 
the temperature and co near the critical point of the medium 
for a specific material characterized by some characteristic 
spectrum. On the other side - when studying critical phe¬ 
nomena one normally starts with some effective Hamiltonian 
where only few basic features of the critical medium are re¬ 
flected. The current article provides a uniform treatment of 
the contributions due both to the van der Waals forces and 
the critical Casimir forces in the framework of an, unavoid¬ 
ably, relatively simple model in which, however, all the cal¬ 
culations have been done on equal footing. The expression 
that we derived for the Hamaker term is, nevertheless, in full 
agreement with the Dzyaloshinskii-Fifshitz-Pitaevskii theory 
E2I33. On the other side the mean-field theory is considered 
as a reliably theoretical horse for the qualitative description of 
the critical phenomena. 
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Appendix A: The Hamaker term for a van der Waals system of 
two different substrates separated by fluid medium 

In this Appendix we derive an expression for the Hamaker 
term for a van der Waals system of two different substrates 
separated by fluid medium. We follow the same lines of action 
when performing the calculations as in Ref. Ea with the 
important generalization that now, within our model (see Sec. 
m the substrates on both sides of the fluid are allowed to be, 
in general, different from each other. 

For the constituent components of the Hamaker constant 
one explicitly derives 

Ai(T,n) = -C(d,o)J l pl(T,ii) < 0, (Ala) 

A Sl ,s 2 {T) = —C(d,<y)J Sl ’ S2 p Sl {T)Ps 2 {T ) < 0, (Alb) 

A Si j{T 7 p) = C(d,o)J Si,l p Sj (T)p b (T,p) > 0, i = l,2, (Ale) 

where C{d , o) = G(d , cr)/((J — 1), and for the direct substrate- 
substrate potential we have assumed, in accordance with Eqs. 

and s that it is given again by van der Waals type 
interaction 


(A2) 

Let us note that A/ < 0 and A 5l ^ < 0, i.e., the direct substrate- 
substrate interaction as well as the fluid-fluid interaction lead 
to attraction between the plates bounding the system while 
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the fluid-substrate interactions are contributing to a repulsion 
between them. Furthermore, let us stress that A/ and A Si i = 
1,2, are singular functions on the temperature near T c being 
dependent on the bulk order parameter Pb , while A < t i2 is an 
analytic function of T near T = T c . 

Using Eqs. ( |A 1 a[ ) - ( |Alc| >, for 11,\ one immediately derives 
from Eq. ( |2.5[ i that 


H f 


’ A (T,p)=-C(d,o){j l pl(T,p)+J^p Sl (T)p S2 (T) 

- Pb {T lP ) [j^ l Psi (T)+J S2 ’ l Ps 2 (T)] } . (A3) 


For the sake of simplicity in the next lines of this section we 
are going to omit giving explicitly the arguments of H A , Pi, 
and p Sj , i = 1,2, as they do not change and are of no relevance 
for the following discussion. Within our model A/ and A s . j, 
i = 1,2 being a part of the grand canonical potential Q. have 
to be proportional to k%T. Obviously, similar statement has to 
hold for A. ?1 . Thus, one has that this shall be also true for 

the pair potentials w XA between any two substances X and 
Y belonging to the system. According to the general theory 
of the dispersion forces between molecules B180II , neglecting 
the quantum effects, for the "zero frequency contribution" to 
w xy (r) one indeed has (see Eq. (6.24) on p. 121 in Ref. 

GMb 


^'" (r) = “2(4^)V a0za0r ’ (A4) 

where Oqx and ctor are the static electronic polarizabilities of 
each of the two substances X and Y. The last implies that 
the constants A/, A i; /, i = 1,2, and ,4 V| Xn are not independent. 
From Eq. ( |A4| > one derives that 

J s '-^J l =J s '-lj s 2- 1 . (A5) 


Thus, it is convenient to express also H A in terms of these 
variables. One obtains that 


/3(ct - 1 )H a = 

_ 4A j—i 

G(d,o)K*x r 


1 

4 


G(d,o)K c Xtyb 


(A6) 


Since Sj = (K/K c )s ix , i = 1,2, with K = fiJ l ST and K, = /3 c J l SI 
it is clear that s b i = 1,2 are temperature dependent, while 
s ijC , i = 1,2 are not. When d = <J = 4 Eq. ( |A6[ > simplifies to 

_ 32 K 

mA = -^k SX - cS2 ' c 

K 7T^ 

+ t^( s i,c + 'S2,c)06 — (A7) 

K c oZ 


where we have taken into account that in mean-filed approxi¬ 
mation p c = 1/2. In Eq. ( |A7| ) the term proportional to A -1 
corresponds to the direct interaction between the confining 
walls, the second (proportional to A 0 ) corresponds to the in¬ 
teraction between the confining walls and the fluid medium, 
while the last one (proportional to A) reflects the interaction 
between the constituents of the fluid medium. The corre¬ 
sponding result for a d = <7 = 3 system is 


A- — ( s i.c4~S2,c) ( l ) b—-z^K(j) b , (A8) 

J\ c O 

which, taking into account Eq. m means that the Hamaker 
constant is equal to 


/^Anam 


n 2 3 4 5 6 7 8 9 XK ]~J 

(= 1,2 


6s, -,c 
nX K, 



(A9) 


Using this equation one can rewrite Eq. ( |A3| ) in the form given 
in Eq. ( |3.16| ) in the main text. There, when calculating the 
critical component of the force, we often characterize the sys¬ 
tem by the values of the parameters A and s,- c , i = 1,2 of the 
parameters s,-, i = 1,2 at the bulk critical point of the fluid 
medium (j3 = j3 c ,Ap = 0) [see Eqs. (3.7 1 , (3.10) and ( |3.11[ )]. 


One must pay attention when using Eqs. (A7 1 and (A8 i, be¬ 
cause the substrate-fluid coupling parameters s, iC , i = 1,2 de¬ 
pend on the dimensionality d of the system and the decay 
exponent a [see Eqs. ( |3. 7 [ > and ( |3. 11 [ >] . The behaviour of 
Eq. ( |A8| > as a function of the temperature and field scaling 
variables x, and je fl is depicted in Fig. [5] 
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